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In a solid material strong interactions between 
the electrons can lead to surprising properties. 
A prime example is the Mott insulator, where 
the suppression of conductivity is a result of in- 
teractions and not the consequence of a filled 
Bloch band^. The proximity to the Mott in- 
sulating phase in fermionic systems is the ori- 
gin for many intriguing phenomena in condensed 
matter physics^, most notably high-temperature 
superconductivity^^. Therefore it is highly desir- 
able to use the novel experimental tools devel- 
oped in atomic physics to access this regime. In- 
deed, the Hubbard model"^, which encompasses 
the essential physics of the Mott insulator, also 
applies to quantum gases trapped in an optical 
lattice^ ^. However, the Mott insulating regime 
has so far been reached only with a gas of bosons^, 
lacking the rich and peculiar nature of fermions. 
Here we report on the formation of a Mott insu- 
lator of a repulsively interacting two-component 
Fermi gas in an optical lattice. It is signalled 
by three features: a drastic suppression of dou- 
bly occupied lattice sites, a strong reduction of 
the compressibility inferred from the response of 
double occupancy to atom number increase, and 
the appearance of a gapped mode in the excita- 
tion spectrum. Direct control of the interaction 
strength allows us to compare the Mott insulating 
and the non-interacting regime without changing 
tunnel-coupling or confinement. Our results pave 
the way for further studies of the Mott insulator, 
including spin ordering and ultimately the ques- 
tion of d-wave superfluidity ^ ^. 

The physics of a Mott insulator is well captured by 
the celebrated Hubbard model which is widely used to 
describe strongly interacting electrons in a solid. It as- 
sumes a single static energy band for the electrons and 
local interactions, i.e. spin- up and spin-down fermions 
are moving on a lattice and interact when occupying the 
same lattice site. The consequence of strong repulsive 
interactions is that even fermions in different spin states 
tend to avoid each other. In the case of a half filled 
band the particles get localised and an incompressible 
state with one fermion per site forms. Since no symme- 
try is broken, the transition between the metallic and the 
Mott insulating regime at finite temperature exhibits a 
crossover rather than a phase transition. 

The Hubbard model ignores various complexities of 
material^ but it has been highly successful to study the 
nature of the Mott insulating regime, including magnetic 



phenomena^ and high-temperature superconductivit}^. 
Yet, despite its simplicity, it turned out that the 
fermionic Hubbard model is in many cases computa- 
tionally intractable and important puzzles remain to be 
solved. In particular, the question whether the ground 
state of the lightly doped 2D Hubbard model supports 
d-wave superconductivity is as yet unanswered. 

Compared to real materials, a fermionic quantum gas 
trapped in an optical latti ce is a much purer realisation 
of the Hubbard modePEI^. It offers a new approach to 
understand the physics of strongly correlated systems. In 
an optical lattice three mutually perpendicular standing 
laser waves create a periodic potential for the atoms. The 
kinetics of the atoms is determined by their tunnelling 
rate between neighbouring lattice sites and the interac- 
tion is due to inter-atomic collisions occurring when two 
atoms are on the same site. In a gas of fermions in differ- 
ent spin states this collisional interaction can be widely 
tuned through a Feshbach resonance without encounter- 
ing significant atom losseJ^^. 

A landmark result has been the observation of the su- 
perfluid to Mott insulator transition using bosonic atoms 
trapped in an optical lattice^. Yet it is the fermionic 
character combined with repulsive interactions which 
provides the intimate link to fundamental questions in 
strongly correlated electron systems. Whereas exper- 
imental studies of fermionic quantum gases in three- 
dimensional optical lattices have been scarce and fo- 
cused on non-interacting and attractively interacting 
cases-^^ we investigate the repulsive Fermi-Hubbard 
model and its paradigm, the Mott insulator. 

In optical lattice experiments the presence of an un- 
derlying harmonic trapping potential has an important 
influence on the observable physics. Let us first con- 
sider a zero temperature Fermi gas prepared in an equal 
mixture of two non-interacting spin components. All 
available single particle quantum states will be filled up 
to the Fermi energy and, for a sufficiently large num- 
ber of trapped atoms, a band insulating region with two 
atoms per site appears in the trap centre, surrounded by 
a metallic shell with decreasing filling, see figure 1. An 
important quantity to characterise the state of the sys- 
tem is the fraction D of atoms residing on lattice sites 
which are occupied by two atoms, one from each compo- 
nent. For the non-interacting case this double occupancy 
should increase in a continuous fashion with the number 
N of atoms in the trap. 

A very different scenario can be anticipated for a gas 
with increasingly st rong repulsive interactions. A Mott 
insulator will appeaif^^CSl, at first in those regions of the 
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FIG. 1: Energy spectrum of a Fermi gas in an opti- 
cal lattice with an underlying confining potential. In 

the non-interacting case (a) the curvature of the lowest Bloch 
band reflects the harmonic confinement. At zero tempera- 
ture all states up to the chemical potential /i are filled with 
atoms of both spin states (green and blue), b, In the Mott 
insulating limit the energy cost for creating doubly occupied 
sites greatly exceeds the temperature T and the kinetic en- 
ergy parametrised by J, giving rise to a gap of order U. The 
energy spectrum of single particle excitations is then depicted 
by two Hubbard bands. Doubly occupied sites correspond to 
atoms in the upper Hubbard band. 



trap where the local filling is approximately one atom 
per site. For very strong repulsion the entire centre of 
the trap contains a Mott insulating phase and double 
occupancy is suppressed, see figure 1. Since the Mott in- 
sulating region is incompressibld^^^, the suppression of 
double occupancy should be robust against a tightening 
of the trapping potential, or equivalently, against an in- 
crease in the number of trapped atoms. However, once 
the chemical potential /i reaches a level where double 
occupation of sites becomes favourable, a metallic phase 
appears in the centre and the double occupancy increases 
accordingly The energy spectrum in the Mott insulat- 
ing phase is gapped, with a finite energy cost required to 
bring two atoms onto the same lattice site. This energy 
has to be large compared to the temperature in order to 
keep the number of thermally excited doubly occupied 
sites small. Thermally excited holes in the centre are 
suppressed by the chemical potential /i. 

Our experiment is performed with a quantum degen- 
erate gas of fermionic ^^K atoms, prepared in a balanced 
mixture of two magnetic sublevels of the F = 9/2 hyper- 
fine manifold {F: total angular momentum). Feshbach 
resonances allow us to tune the 5-wave scattering length 
between a = 240 ± 4ao and 810 ± 40ao as well as to pre- 
pare non-interacting samples. Here ao is the Bohr radius. 
The two-component Fermi gas is subjected to the poten- 
tial of a three-dimensional optical lattice of simple cubic 
symmetry. In terms of the recoil energy En — / {2m\^) 
the lattice potential depth Vb is chosen between 6.5 and 
12Er. Here h is Planck's constant, m the atomic mass 
and A = 1064 nm the wavelength of the lattice beams. 
The system is described by the Hubbard Hamiltonian: 
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The onsite interaction energy is given by U and the 
tunnelling matrix element between nearest neighbours 
{ij) by J. The quotient U/{6J) which characterises the 
ratio between interaction and kinetic energy can be tuned 
from zero to a maximum value of 30. The fermionic cre- 
ation operator for an atom on the lattice site i is given by 
cj^, where a G {t^i} denotes the magnetic sublevel. The 
particle number operator is hi = fii^ + n^^, fiia- = cJ^Qa, 
and Ci is the energy offset experienced by an atom on 
lattice site i due to the harmonic confining potential. 

In order to characterise the state of the Fermi gas in 
the optical lattice we have developed a technique to mea- 
sure the fraction D of atoms residing on doubly occupied 
sites with a precision of down to 1%. The experimental 
procedure is as follows. The depth of the optical lat- 
tice is rapidly increased to 30^^^ to inhibit further tun- 
nelling. In the next step, we shift the energy of atoms on 
doubly occupied sites by approaching a Feshbach reso- 
nance. This enables us to specifically address only atoms 
on doubly occupied sites by using a radio frequency pulse 
to transfer one of the spin components to a third, pre- 
viously unpopulated magnetic sublevel. The fraction of 
transferred atoms is obtained from absorption images and 
allows us to deduce the double occupancy. 

The double occupancy as a function of total atom num- 
ber is plotted in figure 2a, where the non-interacting sit- 
uation is compared to the case of strong repulsive inter- 
actions. The former shows the expected rapid increase 
of double occupancy with atom number. A strikingly 
different behaviour is observed in the strongly repulsive 
regime with U ^ J^T^ ja^ where a Mott insulator is ex- 
pected. The double occupancy is strongly reduced to 
values systematically below 2% for small atom numbers. 
This is direct evidence for the suppression of fluctuations 
in the occupation number and for the localisation of the 
atoms. 

In order to experimentally investigate the compressibil- 
ity on entering the Mott insulating regime we determine 
how the double occupancy changes with increasing atom 
number, i.e. we extract the slope dD/dN from curves 
such as shown in figure 2. This slope is a good mea- 
sure of the compressibility = dn/dja in those regions 
of the cloud where the filling n is near unity or larger, 
since n increases with D. We estimate the filling in the 
trap centre for the non-interacting case from the mea- 
sured double occupanc}/^^. It significantly exceeds one 
atom per site, e. g. (n) = 1.4 for TV = 5 x 10"^, = ^Er 
and a temperature T of 30% of the Fermi temperature 

The slope dD/dN is displayed in figure 3 for a wide 
range of interaction strengths. The data shows that we 
access two regimes: For small onsite interaction energies 
U the slope dD/dN is positive and the system is com- 
pressible. Yet for U/h > 5 kHz the measured compress- 
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FIG. 2: Double occupancy in the non- interacting and 
Mott insulating regime, a, A significant increase of tfie 
double occupancy with atom number is observed in the non- 
interacting regime (empty circles) whereas on entering the 
Mott insulating regime the double occupancy is suppressed 
(filled circles) . The corresponding onsite interaction strengths 
are U/h = dz 80Hz and U/h = 5.0 =b 0.6 kHz. b, In the 
Mott insulating regime the double occupancy is strongly sup- 
pressed. It starts to increase for large atom numbers indicat- 
ing the formation of a metallic region in the trap centre. The 
blue and red lines represent the theoretical expectation for D 
in the atomic limit (see text and methods summary). Values 
and error bars are the mean and s. d. of 4 to 8 identical 
measurements. The systematic relative errors for the atom 
number, double occupancy and lattice depth are estimated 
to be 20%, 10%, and 10% respectively, with corresponding 
relative errors in J of up to 30%. These systematic errors 
apply to all further measurements. 



ibility vanishes. This indicates that we have entered the 
Mott insulating regime. It implies a large central region 
with a filling reduced to one atom per site, surrounded 
by a metallic region with lower filling. 

Further insight is gained by comparing our measure- 
ments with the theoretical values of dD/dN calculated 
in the atomic limit^^ of the Hubbard model, including 
confinement and finite temperature. In this limit the ki- 
netic energy is neglected by setting the tunnelling matrix 
element J to zero (see also methods summary). We find 
good agreement between theory (black line in figure 3) 
and experimental data for U ^ 6J, where the above 
assumption is acceptable. For the calculation we have 
assumed a temperature of T = 0.28Ti7, which is deduced 
from the entropy in the dipole trap as discussed in the 
methods section. For zero temperature the slope dD/dN 
would vanish as soon as U becomes larger than the chem- 
ical potential /i, which is h x 2.7 kHz for = 8 x 10^ 
atoms and a lattice potential of Vb = l2Eii. Both our 
measurements and the model at finite temperature show 
a finite compressibility extending beyond U/h = 2.7 kHz, 
which can be attributed to thermal excitations. For the 
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FIG. 3: The transition to an incompressible sample. 

Upon changing U, two regimes can be distinguished by the 
slope DD/dN. For vanishing interaction the large initial slope 
signals the filling of the Bloch band. Increasing U reduces 
the double occupancy. For U/h > 5 kHz a change in atom 
number can no longer change the double occupancy. The 
compressibility dD / dN is obtained from a least squares fit of 
D{N) = {dD/dN)N + Do to data such as shown in figure 2, 
with atom numbers in the interval from 25 x 10"^ to 8 x 10^. 
Error bars denote the confidence interval of the fit. The ex- 
pected slope in the atomic limit is illustrated with a black line 
for a lattice depth of 12Er and T = 0.28Tf. 



largest attained interaction U/h = 8.1kHz the thermal 
excitations are characterised by T = O.llU/k = 0.28Tp 
corresponding to 3% vacancies in the trap centre accord- 
ing to the theoretical analysis presented in the methods 
section {k is Boltzmann's constant). The vanishing slope 
dD/dN at this filling implies incompressibility of the 
core. The obtained ratio T/U is comparable to estimates 
for the bosonic Mott insulatoi^^. 

In the strongly repulsive regime, the measured com- 
pressibility should vanish if ja < U. For atom numbers 
corresponding to higher chemical potentials a metallic 
phase will appear in the trap centre and the double oc- 
cupancy will increase. We observe this characteristic 
behavioui!^ which is a consequence of the presence of 
a Mott insulator, see figure 2b. The behaviour agrees 
well with that expected from the Hubbard model in the 
atomic limit (lines in figure 2b). The free parameters in 
the theory curves, the temperature and a constant offset 
in D^ are determined by a least squares fit to the data. 
The fits yield temperatures of 0.2 ±0.1Tp. However, the 
accuracy is limited due to the high sensitivity to the en- 
ergy gap and the harmonic confinement. The constant 
offset in D accounts for the finite double occupancy in 
the ground state caused by second order tunnelling pro- 
cesses as well as a systematic offset of 0.5% stemming 
from technical imperfections in the initial preparation of 
the spin mixture. 

An important feature of a Mott insulator is the energy 
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FIG. 4: Emergence of a gapped mode. With increasing 
interaction (blue, green, red) the response to modulation of 
the lattice depth shows the appearance of a gapped mode. 
The weight of this peak grows with U / (6 J) and saturates. All 
modulation spectra were obtained with (32 ±2) x 10^ atoms. 
The weight of the peak shown in the bottom right panel is 
X;^Azy [Diui) - ^{D{200}lz) + D{700llz))], where D{iyi) is 
the measured double occupancy at frequencies Ui which are 
evenly spaced in steps of Au = 500 Hz. It is plotted in units 
of 12J/h. The first four datapoints are taken for a lattice 
depth of 6.5Er, the next at 7, 8 and IOEr, from left to right, 
respectively. The lines serve as a guide to the eye. Values 
and error bars are the mean and s. d. of 4 to 8 identical 
measurements. 

gap in the excitation spectrum. The lowest lying exci- 
tations are particle-hole excitations centred at an energy 
U. The actual gap in the energy spectrum is reduced 
with respect to this value due to the width of the en- 
ergy bands experienced by particles and hole^. A suit- 
able technique for probing this excitation spectrum is to 
measure the response of the quantum gas in the opti- 
cal lattice to a modulation of the lattice deptlP^H^ we 
apply 50 cycles of sinusoidal intensity modulation of all 
three lattice beams with an amplitude of 10%. The re- 
sponse is quantified by recording the double occupancy 
as a function of modulation frequency. With increasing 
interactions we observe the emergence of a gapped mode 
in the excitation spectrum (figure [4| . For small values of 
U I (6 J), the double occupancy is not affected by the mod- 
ulation of the lattice depth but for large values of [// (6 J) 
a distinct peak appears for modulation frequencies v near 
U jh. Furthermore, the area under the excitation curve 
divided by 12 J jh as a measure for its width increases 
with interaction strength and starts to saturate at large 



values oil] I (6 J), see figure [4]i. 

The presented approach to the physics of the repulsive 
Fermi-Hubbard model is completely different and com- 
plementary to that encountered in solid-state systems, 
and provides a new avenue to one of the predominant 
concepts in condensed matter physics. In this first ex- 
periment we have found clear evidence for the formation 
of a Mott insulator of fermionic atoms in the optical lat- 
tice. We could set limits for the deviation from unity 
filling in the Mott insulator by directly measuring the 
residual double occupancy and by deducing the number 
of holes from a realistic estimate of the temperature. The 
temperature is found to be small compared to the onsite 
interaction energy and the Fermi temperature. In ad- 
dition, we have obtained good quantitative agreement 
with the Hubbard model in the atomic limit for a wide 
range of parameters. In further investigations of e.g. the 
energy spectra, the high resolution achieved may give 
direct insights into the width of Hubbard bands , the 
lifetime of ex citations and the level of anti-ferromagnetic 
orderin^i^l^ in the system. 



Methods summary 

In the atomic limit U ^ 6J of the Hubbard model 
we assume full localisation of the fermions and thus ne- 
glect the kinetic energy. Each site is treated in the grand 
canonical ensemble with three possible occupation num- 
bers n e {0,1,2}. The partition function for site i is 
then Zi = ^^z^'eM-f^Ei.n) = 1 + 2zexp(-/3eO + 

exp(— 2/3ei — PU) where /3 = 1/kT is the inverse tem- 
perature, z = exp(/3/i) the fugacity, /i the chemical po- 
tential, Ei^n the energy of n particles on site i, and 
the energy offset due to the harmonic confinement. For 
the probability to find a double occupancy {d)i or a va- 
cancy {v)i one obtains {d)i = z^exp(— 2/3e^ — /3U)/Zi 
and {v)i = l/Zi. Double occupancy D and total parti- 
cle number N of the system are obtained by summing 
over all sites, e.g. D = 2{d)i/N^ where the equation 
for N is first solved numerically with respect to z. The 
entropy is = ^(/cT InZ^). We calculate the tem- 
perature in the lattice by assuming that this entropy is 
the same as the entropy determined from temperature 
measurements in the dipole trap (see methods) . The fits 
in figure 2b involve U as determined by modulation spec- 
troscopy {U/h = 4.7 ± 0.1 kHz and 6.1 ± 0.1 kHz) since 
band structure calculations disagree with the measured 
value by up to 30% for the largest scattering lengths. 
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Methods 

Preparation. After sympathetic coohng with ^'^Rb, 
2 X 10^ fermionic ^^K atoms are transferred into a dipole 
trap operating at a wavelength of 826 nm. Initiahy, a 
balanced spin mixture of atoms in the Ittif) = |— 9/2) 
and |— 7/2) states is prepared and evaporatively cooled 
at a magnetic bias field of 203.06 Gauss. Using this mix- 
ture we realise non-interacting samples with a scatter- 
ing length of a = ± lOao. Repulsive interactions are 
obtained by transferring the atoms in the |— 7/2) state 
to the |— 5/2) state during the evaporation, thus cool- 
ing and preparing a spin mixture of atoms in |— 9/2) 
and |— 5/2) states, close to a Feshbach resonance at 
224.21 GausP^. After tuning the scattering length to 
the desired value we load the atoms into the lowest Bloch 
band of the optical lattice by increasing the intensity of 
three retroreflected laser beams within 200 ms using a 
spline ramp. The beams have circular profiles with 1/e^ 
radii of (160, 180, 160) /im at the position of the atoms. 
For a given scattering length and lattice depth J and 
U are inferred from the Wannier functions including the 
interaction induced coupling to the second Bloch band. 
The latter leads to corrections of up to 15% with respect 
to the single band model. 

Radio -frequency spectroscopy. By increasing the 
depth of the optical lattice to 30 Eji in 0.5 ms tunnelling is 
suppressed. The magnetic field is tuned to 201.28 Gauss, 
where a molecular state for a |— 9/2), |— 7/2) pair with 
binding energy x 99 ± 1 kHz and a weakly interacting 
state for a |— 9/2), |— 5/2) pair exist^^. A radio-frequency 
TT-pulse dissociates (associates) pairs and changes the 
spin state of those |— 7/2) (|— 5/2)) atoms that share a 
site with a |— 9/2) atom. Finally the magnetic field is in- 
creased to 202.80 Gauss dissociating any molecules and 
the lattice potential is ramped down in 10 ms. All con- 
fining potentials are switched off and the homogeneous 
magnetic bias field is replaced by a magnetic gradient 
field in the same direction applied for 2 ms, thus spatially 
separating the spin states. 

Imaging. After 6 ms of time-of- flight all three clouds 
are imaged simultaneously. Due to a reproducible change 
of the imaging beam profile between the atomic absorp- 
tion image and the subsequent reference image without 
atoms, residual structures are present in the density pro- 



files. These are reduced by repeating the entire exper- 
iment without loading atoms and subtracting the ob- 
tained residual density distribution from the atomic den- 
sity distribution. The number of atoms Nm^ per spin 
component mi? is determined from the 2D column den- 
sities by simultaneously fitting the sum of three quar- 
tic terms A ■ max(l — {x/wx)'^,0) • max(l — {y /wy)'^,0) 
with identical widths Wx^y and mutual distances. This 
permits accurate detection of atom numbers down to 
200 atoms per spin state. We have validated the ab- 
solute accuracy of the fits against integration of the den- 
sity. The fraction D of atoms residing on doubly oc- 
cupied sites is defined as I) = 2A^^^//7V where N = 
N_g/2 + N_Y/2 + and m^^ = -5/2 (-7/2) for 

samples initially containing atoms in the |— 7/2) (|— 5/2)) 
states, respectively. The relative uncertainty in D is 10%, 
validated against m easur ements of the adiabatic molecule 
formation efficiency^^^^l We estimate the relative sys- 
tematic error for the total atom number N to be less 
than 20%. The |-9/2), |-5/2) mixture shows an offset 
of 0.5% in D due to |— 7/2) atoms remaining from the 
initial spin transfer during evaporation. 

Temperature. The temperature is measured in the 
harmonic dipole trap before ramping up the lattice and 
after a subsequent reversed ramp back into the dipole 
trap. The highest temperatures measured before and af- 
ter ramping are = 0.1 5Ti? and Tf = 0.24T;?, respec- 
tively. Since we expect non-adiabatic heating to occur 
during the lattice ramp up as well as during ramp down, 
we use the mean value of 0.195Tp as a realistic estimate. 
With this, we calculate a temperature of T = 0.28Tp 
in the Mott insulating regime (a = 810ao, Vq = 12Eji, 
N = 10^), corresponding to 3.3% holes and a compress- 
ibility as low as dn/dfi = 0.09//i in the centre. For the 
temperatures in the dipole trap before and after the lat- 
tice ramp we would obtain 0.3% holes for Ti and 11.5% 
for Tf. The reported temperatures represent upper lim- 
its, since we have achieved temperatures down to 0.08Tp 
in the dipole trap prior to loading. Due to inelastic colli- 
sions we lose at most (4.8 ±0.6)% of the atoms during the 
preparation of the Mott insulating state for the parame- 
ters above, where the losses are expected to be highest. 
The inelastic decay time for atoms on doubly occupied 
sites exceeds 850 ms, which is significantly longer than 
the relevant experimental timescale. 



